Charge conservation and the shape of the ridge of two-particle correlations in 

relativistic heavy-ion collisions 
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We demonstrate that in the framework of the event-by-event hydrodynamics followed by statistical 
hadronization, the proper charge conservation in the mechanism of hadron production provides the 
crucial non-flow component and leads to agreement with the two-dimensional two-particle correlation 
data in relative azimuthal angle and pseudorapidity at soft transverse momenta [pr < 2 GeV). The 
fall-off of the same-side ridge in relative pseudorapidity follows from the fact that a pair of particles 
with balanced charges is emitted from the same fluid element, whose collective velocity coUimates the 
momenta of the pair. We reproduce basic experimental features of the two-dimensional correlation 
function, such as the dependence on the relative charge and centrality, as well as the related charge 
balance functions and the harmonic flow coefficients as functions of the relative pseudorapidity. 
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Two-particle correlation functions in the relative angle 
A</) and pseudorapidity Ary are valuable tools to study 
collective flow and the mechanism of particle emission 
in relativistic heavy-ion collisions. The harmonic com- 
ponents of the collective flow are visible in the dihadron 
correlation function as two ridge structures on the same 
(A</) ~ 0) and away (A(^ ~ tt) sides [l|]. There is, how- 
ever, an on-going discussion concerning the puzzling na- 
ture of the same-side ridge While it is commonly 
accepted that the collective harmonic flow [l| determines 
the profile in A(/) for large pseudorapidity separations, up 
to now the shape of these structures in A77, in particular 
the rather fast fall-off of the same-side ridge, remains an 
object of active debate, with arguments that the presence 
of (mini)jets Q is essential to explain the phenomenon 
and that the applicability of hydrodynamics, reproduc- 
ing numerous other features of the heavy-ion data, is at 
stake. Thus the issue is of great importance for the funda- 
mental understanding of relativistic heavy-ion collisions. 
Other attempts to explain the nature of the ridge can be 
found in Refs. Q. 

In this Letter we show that two basic features of 
the two-particle correlations get a quantitative explana- 
tion via the charge balance mechanism of particle emis- 
sion: 1) the shape of the same-side ridge in A77, and 
2) the difference between the correlation functions for 
like- and unlike-sign particles. Thus we explain the 
ridge puzzle in a natural way, amending the (event- 
by-event, 3-|-l-dimensional, viscous) hydrodynamics with 
the local charge- conservation mechanism in the statisti- 
cal hadronization occurring after the hydrodynamic evo- 
lution. This important charge balancing [5rt8j, simply 
stating that the hadron production conserves locally the 
charge, is an otherwise well-known and measured feature. 

The results presented in this work concern "soft 



physics" (typically with the transverse momentum of all 
particles pt < 2 GeV) and unbiased correlations, where 
the kinematic cut on both particles is the same. The 
relevant correlation function is determined as 



C(Ar?,A0) 



<:i(A^,A</,)/7VP-,(A,y,A0), (1) 



where ^reaTmixcd(A'?i A(/)) denote the two-dimensional 
distributions of pairs of particles with relative pseudo- 
rapidity A77 and azimuth Ac/), obtained from the real 
and mixed events, respectively. Our approach consists 
of using GLISSANDO [9] to generate the Glauber-model 
initial condition, then running event-by-event 3-l-lD hy- 
drodynamics with shear and bulk viscosities 10], and 
finally carrying_^out the statistical hadronization with 
THERMINATOR [llj at the freeze-out temperature T/. Our 
simulations incorporate the kinematic cuts of the STAR 
experiment, with I77I < 1, appropriate px cuts specified 
later, as well as the detector efficiency at the level of 




FIG. 1. (Color online) A schematic view of the charge bal- 
ancing mechanism, producing pairs of particles with opposite 
charges. The rectangles indicate fluid elements moving out- 
ward with a collective velocity u. The dot indicates the space- 
time location of the emission of the pair of opposite-charge 
particles of momenta pi and P2- The dashed line represents 
a neutral resonance, decaying into a pair particles. 
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FIG. 2. (Color online) The correlation function C for like-sign (a,c) and opposite-sign (b,d) pairs. Panels (a,b) and (c,d) 
correspond to absent and present direct charge balancing, respectively. Inclusion of charge balancing sharpens the peak around 
A77 = = and causes the desired fall-off of the same-side ridge (centrality 30 - 40%, Tf = 140 MeV, 0.2 < < 2 GeV). 



90%, estimated to hold for the registered charged parti- 
cles in STAR. For simplicity, we set all chemical poten- 
tials at freeze-out to zero, which is a good approximation 
at RHIC. Other, more technical details of our approach 
may be found in Ref. [To|). 

The observed charge balance functions can be ex- 
plained assuming that opposite charge pairs are created 
towards the end of the evolution 0, To implement 
this mechanism in a simple model way but with a real- 
istic hydrodynamic flow (we call it direct charge balanc- 
ing), we enforce that the same-species charged hadron- 
antihadron pairs are produced at the same space-time 
location x (see Fig. [T]) . The hadron momenta pi and p2 
are determined independently according to the Cooper- 
Frye formula. The fact that the fluid element moves with 
a collective velocity u^^[x) causes a certain-degree of col- 
limation of the momenta of the produced pair. An ad- 
ditional balancing mechanism comes from the decays of 
neutral resonances (see Fig.[T]). The correlations induced 
by balancing are of a non-flow character, i.e., cannot be 
obtained by the folding of single-particle distributions 
containing the collective flow. 

To illustrate the relevance of the effect, in Fig. [5] we 
show the results of our simulations for several cases for 

the like-sign (H — h, ) and unlike-sign (H — ) pairs. In 

panel (a) we show the correlation C(-|--|-, ) without 

direct balancing. We note the completely flat ridges, re- 
flecting the approximate boost-invariance in the inves- 
tigated kinematic range and, of course, the presence of 
flow. We use the framework of event- by- event viscous hy- 
drodynamics which generates realistic elliptic and trian- 
gular flows in the collisions 12]. Therefore the dominant 



modulation of the shape in azimuth of the elliptic and tri- 
angular flows is well reproduced [l^ . Panel (b) shows 
the same for C(-| — ), where some mild fall-off in A?] of 
the same-side ridge follows from the resonance decays. 
Panels (c) and (d) include the direct charge balancing. 
We now note a prominent fall-off of the same-side ridge 
in C{-\ — ), which is our key observation: the quantity 
C(H — ) — 1 drops from the central region to jAjyl = 2 
by about a factor of 2. The fall-off is also enhanced 

for C(++, ) due to secondary effects from balancing 

of heavier particles, which later decay. The results of 
panels (c,d) are in qualitative and approximate quanti- 
tative agreement to the results of the STAR Collabora- 
tion, where the HBT correlations for identical particles 
are subtracted j2j . 

To check whether our mechanism is correct also at 
the quantitative level, we now proceed to the inves- 
tigation of the charge balance functions, defined as 
B(A77) = ( - N++)I{N+) + {N_+- N__)/{N_), 
where {Nab) denotes the event- aver aged distributions of 
particles a and h with relative rapidity Ary, and {Ng) 
stands for the average number of particles a in the 
acceptance window It^I < 1. We note that the charge 
balance function is related to the distributions in the 
numerator of Eq. ([l} , 



B(Ar/) 



/ dA0[Arfi''(A?7, A</>) - iV57(A77, A(/.)] 



(2) 



The outcome, with correct agreement to the data, is pre- 
sented in Fig.[3l We note a preference to the lower freeze- 
out temperature, Tf = 140 MeV. 
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FIG. 4. (Color online) The flow coeflicients w| (top lines) and (bottom lines). Simulations with direct charge balancing 
are drawn with thick solid (T/ — 140 MeV) and dashed lines (T/ = 150 MeV), while the corresponding reference simulations 
without direct charge balancing are drawn with thin lines. The dashed bands are extracted from the fits to experimental data 
reported by STAR in Table I of Ref. (0.15 < pt < 4 GeV, as in the experiment). 



The next quantitative investigation concerns the de- 
pendence of the flow coefRcients on Ar/, defined as 

v2(A?7)= J dA(j)cos{nA(j))C{Ar],A(j)). (3) 

The projection on the A77 axis of the different harmonics 
yields the squares of the consecutive flow components u„ 
present in the dihadron correlation functions. The results 
presented in Fig.|3]show agreement with the experiment, 
best for the mid-peripheral collisions and Tf = 140 MeV. 
For the peripheral collisions, where the hydrodynamic ap- 
proach is less justified, the agreement is qualitative, indi- 
cating that the hydrodynamic calculation overestimates 
the elliptic flow for large centralities. The experimental 
bands are extracted by integrating a model function fit 
to the measured dihadron correlations j2!| , varying the fit 
parameters within the estimated uncertainty. We note 
that these uncertainties are large for the central and pe- 
ripheral cases. Our simulations incorporating the direct 
charge balancing (thick lines) exhibit the quested fall-off 
with |A?7|, while the cases without direct balancing (thin 
lines) are flat. The independence of on A77 for the 
emission without charge balancing reflects the approxi- 
mate pseudorapidity independence of the collective flow 
in the considered kinematic window. Charge balancing 
induces an additional component in C, of limited range 
A77I ~ 1. The coUimation of the opposite charge pairs 
occurs in the relative angle as well 

BUI- As a result. 



the contribution from charge balancing in C{Ai], A0) ac- 
quires the form of a 2-dimensional peak at Ar] = A4> = 0. 
The shape in Ar/ of the non-flow component in w| is 
qualitatively reproduced in the simulations, but the over- 
all strength is somewhat larger than extracted from the 
model fit in [2]. Thus our study shows that the charge 
balancing is the non-flow source of the observed A77 de- 
pendence of the flow coefficients JLSj. The qualitatively 
similar behavior of higher-order harmonics, which needs 
higher statistics in our simulation, as well as w^, where 
the effects of the transverse- momentum conservation (not 
included in the present study) are important [Tcj , will be 
presented elsewhere. 

One may also compare the correlation function C di- 
rectly to the data shown, e.g., in Figs. 1 and 2 of Ref. [l^l , 
obtained for 0.8 < < 4 GeV, and with the HBT 
peak for the same-sign pairs removed. Our simulations 
shown in Fig. [S] display, for the first time in an approach 
based on hydrodynamics, all qualitative features of the 
data and remain also in fair quantitative agreement. In 
particular, we note the proper dependence on the rel- 
ative charge and centrality. Notably, the combinations 

C(+— ) — C(++, ) obtainable from Fig. [S] exhibit no 

ridges whatsoever, as they cancel out, leaving the central 
peak as the only structure. 

In conclusion, we remark that the presented sim- 
ple effect based on the local charge conservation in 
the hadronization process is generic in its nature. It 
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FIG. 5. (Color online) Our simulations for the correlation function C with direct charge balancing included for the like-sign 
(a,b,c) and unlike-sign (d,e,f) pairs at three sample centralities (T/ = 140 MeV, 0.8 < pt < 4 GeV as in Ref. (17j]). 



should manifest itself in all approaches where the charge 
balancing is combined with a collective motion of the 
source. Our approach, based on the fluctuating Glauber- 
model initial conditions, state-of-the-art hydrodynam- 
ics, and statistical hadronization incorporating the direct 
charge balancing, is capable of reproducing all basic fea- 
tures of the data for the unbiased correlation function 
C(A7/, A0), as well as for the related quantities, such as 
the charge balance function and the harmonic flow coef- 
ficients w^(A77). The correlation from charge balancing, 
yielding a two-dimensional central peak, comes on top of 
the ridge structures following from the presence of the 
azimuthally asymmetric collective flow It thereby 



brings in a crucial non-flow component in the harmonic 
flow coefficients w^, with a characteristic fall-off in the 
relative pseudorapidity. Thus the collective flow together 
with the local charge conservation is the key to a suc- 
cessful explanation of the shape of the correlation data 
in relativistic heavy-ion collisions. 

We thank Paul Sorensen for a discussion concern- 
ing the fall-off of the ridge, and to Mikolaj Choj- 
nacki for help in modifying THERMINATOR. Supported 
by Polish Ministry of Science and Higher Education, 
grant N N202 263438, and National Science Centre, 
grant DEC-2011/01/D/ST2/00772. The simulations 
were partly carried out on the Cracow Cloud One cluster. 



[1] J. Takahashi, B. Tavares, W. Qian, R. Andrade, 

F. Grassi, et ai, .Phy s.Rev.Lett. 10 3 242301 (2009), 
M. Luzum, [Pi^^s.Lett. B6 96, 499 (201lj] 

[2] G. Agakishiev et al. (STAR), (2011), 

]axXl^:1109.4380 [nucTe^ _ 
[3] T. A. Trainor, (2011), arXiv:1109.2540 [hep-ph] 
[4] S. Gavin, L. McLerran, and G. Moschelli, 

|Phys. Rev. C79 , 051902 (200 9)^ S. Gavin and 

G. Moschelli, |Phys. Rev. C85, 014905 ( 2012) 

[5] S. A. Bass, P. Danielewicz, and S. Pratt, 

Phys. Rev. Lett. 85,1:689"(2000) 
[6] S. Jeon and S. Pratt, Phys. Rev. C65, 044902 (2002) 
[7] P. Bozek, W. Broniowski, and W. Florkowski, Acta 

Phys. Hung. A22, 149 (2005). 
[8] M. M. Aggarwal et al. (STAR), 

|Phys. Rev. C82, 024905 (2010)' 
[9] W. Broniowski, M. Rybczyhski, and P. Bozek, 

Comput. Phys. Commun. 180, 69 (2009) 
[10] P. Bozek, Phys.Rev. C85, 014911 (2012); P. Bozek and 



W. Broniowski, Phys. Rev. C85, 044910 (2012). 

[11] M. Choj nacki, A. Kisiel, W. Florkowski, and W . Bro- 
niowski, I Comput.Phys.Commun. 183, 746 (2012) [. 

[12] B. Schenke, S. Jeon, and C. Gale, 

J.Phys.G G38, 124169 (201l)j Z. Qiu, C. Shen, 
and U. W. Heinz, Phys.Lett. B707, 151X2012); B. H. 
Alver, C. Gombeaud, M. Luzum, and J.-Y. OUitrault, 
|Phys. Rev. C82, 034913 (2010) | 

[13] K. Werner, I. Karpenko, M. Bleicher, 
T. Pier og, and S. Porteboeuf-Houssais, (2012), 
arXiv: 1203.5704 [nucl-th] 

[14] P. Bozek, Phys.Lett. B60 9, 247 (2005) 

[15] K. Aamodt et al. (ALICE Collaboration), Phys.Rev.Lett. 
107, 032301 (2011). 

[16] N. Bo rghini, P. M. Dinh, and J.-Y. OUitrault, 
[Phys. Rev. C62, 034902 (2000); A. Bzdak, V. Koch, 
and J. Liao, Phys.Rev. C83, 014905 (2011) 

[17] B. Abelev et al. (STAR Collaboration), (2008) 
,arXiv:0806.0513 [nucl-ex], 



